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Abstract 

We study Cqo with the use of Thomas-Fermi theory. A spherical shell model 

is invoked to treat the nuclear potential, where the nuclear and core charges 

are smeared out into a shell of constant surface charge density. The valence 

electron distribution and the electrostatic potential are efficiently computed 

by integration of the Thomas-Fermi equation, subject to the shell boundary 

conditions. Total energy is numerically calculated over a range of shell radii, 

and the mechanical stability of the model is explored, with attention to the 

constraints of Teller's theorem. The calculated equilibrium radius of the shell 

is in good agreement with experiment. 
PACS Numbers: 31.15.Bs, 36.40.Qv, 31.10.-Fz 
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The highly symmetrical structure of Cgo has motivated geometrical approximations which 
have previously been invoked to study electronic and optical properties of the molecule. 
While consideration of the icosahedral structure of the molecule is necessary for detailed 
comparison with experiment, previous studies have had success in describing some of the 
properties of Cgo within the continuum approximation, where a system of free electrons are 
constrained to moving on the surface of a sphere 

While the peak electron density should be found on the shell, electrostatic consideration 
of the mechanical stability of the entire system requires that a sizeable fraction of the total 
number of valence electrons be inside the shell. Motivated by this observation, we study 
here a generalization of the previously considered continuum model: we allow the valence 
electrons to move in three dimensions in the external potential generated by a spherical 
shell of constant surface charge density. We call such an artificial molecule whose nuclear 
potential has spherical shell symmetry "spherene." 

We treat spherene by the Thomas-Fermi (TF) method. While TF results are typically 
rather rough, it is often used to efficiently generate starting potentials for more exact self- 
consistent field methods. We have had success [Q using the resulting TF potential in this 
fashion. Of course, TF theory has historically had value in its own right. In this Letter, 
we use the TF results to discuss the stability of Ceo- This is rather subtle business, as it is 
well-known that local density methods such as TF are cases for Teller's theorem 0] which 
states that molecules in TF theory will unbind. We prove Teller's theorem in the context of 
a spherical shell, and we circumvent it by considering the true point charge distribution in 
the calculation of the "nuclear" energy. 

The approach used here was previously employed by N. March to investigate molecules 
with the form of a central atom with tetrahedrally or octahedrally coordinated ligands, such 
as in the case of CH4 and SFg. We closely follow March's work, applying it to the case of Ceo, 
where there is no central atom and an icosahedral arrangement of "ligands." The absence 
of a central atom changes the boundary condition of March at the origin. The case of Ceo 
would seem to be ideally suited to this approach, given its high "coordination" number; and 
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the case of an endohedrally-doped fullerene, such as La@C6o, is of the exact form considered 
by March. 

We start with a positively charged spherical shell of radius R and charge Ze. The shell 
charge arises from the sum of the positive nuclear charges with the core electrons of the 
constituent atoms. We take for the case of Ceo, Z = 60; thus, the valence electrons are the 
remaining 60 vr-electrons. These valence electrons interact with the shell via a spherically 
symmetric cut-off Coulomb potential, given by 



Ze/R 0<r <R 

(1) 

Ze/r r > R 

We can view this model as arising from an approximation of the true nuclear potential 
where all but the monopole term is neglected. The validity of the spherical shell model 
can be examined by expanding the nuclear potential in multipole moments. We denote the 
location of the i*^ atom by a radius R and a set of spherical angles Qi = {9i,(j)i). We center 
our coordinate system on the geometric center, and we align our axes with the five-fold and 
two-fold axes of the molecule. Thus, for the region external to the cage. 
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Ur) = E ^J^^Q^^YLm (2) 



Qim is the 2^-pole moment, given by 



Qim = E Ye^m (3) 

The summation in Eq. § is only over even £ as a result of the inversion symmetry of Ceo- 
Furthermore, in general we note that Qem is non- vanishing only if the spherical irreducible 
representation, denoted by £, when decomposed in terms of the irreducible representations 
of Ih contains the trivial (aig) representation. 

Thus, after the non- vanishing monopole moment, the next non- vanishing elements are 
in £ = 6, followed by £ = 10. For Ceo, we need not consider £ > 10, as the highest lying 
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electron orbital is derived from an £ = 5 manifold. Qim is also only non-vanishing for m = 
and ±5. 

We estimate the error of neglecting ^ 7^ terms by evaluating the relevant dimensionless 
parameters 



47r Q 



(4) 



2e + 1 R'Qoo 

We find that ae.o = 0.026, ae.s = 0.020, aio,o = 0.021, and 010,5 = 0.034. As aim < 1 for 
i < 10, we conclude that the spherical approximation is reasonable for our purposes. 

We consider the dimensionless TF equation without exchange effects at temperature 
T = 0, 



1 

X2 



(5) 



X is the distance from the center of the shell in units of 



6=i 
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where Oq is the Bohr radius of hydrogen. is related to the potential in the usual way 

Ze 



V{r) 



-(f){x) 



(7) 



Without nuclear charge at the origin, the standard atomic boundary condition at x = 
is altered to 0(0) = 0, as the potential is now finite at the origin. The presence of the shell 
gives rise to a discontinuity in the derivative of at the shell. Thus, 



1 

X 



(8) 



where X is the shell radius in dimensionless units and differentiation is with respect to x. 
Additionally, itself is continuous over its domain, and — > as a; -h- 00. 

We obtain numerical solutions to Eq. ^ subject to the above boundary conditions for 
different values of X. A variation of the shooting method |^ is used where we choose a 
trial slope for at the origin, and we integrate outward to the asymptotic region (x ^ X). 



The boundary condition at infinity is replaced by requiring that vanish at an outer shell 
of large radius. The slope is subsequently varied in a systematic way until the boundary 
condition on the outer shell is satisfied. 

The following identities are used as a final check of the numerical procedures: (1) con- 
servation of electron number requires that satisfy 
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(f)2x2dx = 1 (9) 
and (2) the virial theorem for shell systems requires that 

(10) 



9 /"^ , 5 _1 , f°° 02X2 

)2x 2rfx+ / dx = —2 



X 
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where x> is the larger of x and X. The resulting solutions obey the above relations to an 
accuracy of better than three parts in 10"^. 

In Fig. we show the electron charge density n as a function of x obtained from our 
solution for for parameters corresponding to those modeling Cgo, where X = 29.7592 
[R = 6.73ao) and Z = 60. While the potential is found from Eq. n is computed from the 
relation 
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We note that (and consequently n and V) is strongly peaked at the shell, in a consistent 
fashion with the continuum models which constrain the valence electrons to the surface of 
the shell; however, it is significant that nearly 43% of the valence electrons are contained 
inside the shell. We return to this point in our discussion of stability. 

We follow March [0 and conclude that the electronic energy Eg can be simplified to a 
form requiring only values of and its derivative evaluated just inside the shell, 

= f^t^*^^^"^ ~ - 3] (12) 

For the Ceo parameters, we find Ee = —486 Ry. 

Using the Hellman-Feynman theorem, we calculate the radial force that the electrons 
exert on the shell, F^: 



dR 

= -z/-^^ (14) 

dr 

= -^,[X'P\X-)-<P{X-)] (15) 

Only electrons in the interior of the shell can exert a force on the shell, a consequence of 
Gauss' Law. Hence, the presence of charge in the interior provides a centripetal, stabilizing 
force which opposes the centrifugal self-force of the shell. 

From dimensional considerations, we write the self-interaction energy of the shell as 

Un = (16) 

For the uniform shell, c = |. Within TF theory. Teller's theorem implies that the 
stabilizing force of the electrons on the shell is of insufficient magnitude to compensate for 
the repulsive self-force of the shell. Thus, there is no finite equilibrium radius for the shell for 
c = |. We sketch a proof specific to the shell. The proof proceeds by reductio ad absurdum. 
We assume that there exists a finite equilibrium shell radius Rq. Using c = |, Eqs. ^ 



and Eq. 16 imply that 



E{Ro) = ^ZeV{R^) (17) 

Since V is positive for finite r, we conclude that E{Ro) is positive as well. However, the 
virial theorem states that at equilibrium 

E{Ro) = -T (18) 

where T is the total kinetic energy of the electrons. Since T must be positive, we conclude 
from Eq. |18| that E{Rq) is negative. But this is in contradiction to the result of Eq. p!7| . 
Hence, there is no finite R^. 

The continuum approximation overestimates the shell self-force. If one computes the 
self-force by considering the point structure of the ions, a stable equilibrium is obtained. 
For the system of Z ions of charge e located on the vertices of a truncated icosahedron of 
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radius R, we find c ~ 0.4311. We use this value of c to compute the total energy of the 
system, E = Ee + Un, at different shell radii. The resulting energy curve is plotted in Fig. |^. 
From the minimum of the curve, we extract an equilibrium radius, Rq = 7.36ao, which is in 
good agreement with the experimental value of 6.73ao. 

It is amusing to observe that at equilibrium, c is precisely the fraction of valence electrons 
contained inside the shell. This condition follows from simple electrostatic considerations. 

In addition to giving an equilibrium radius in good agreement with experiment, TF 
gives a potential which is an excellent starting potential for more rigorous self-consistent 
field techniques. Furthermore, there are many enhancements of this method which can be 
easily incorporated, such as the inclusion of exchange and correlation and density gradient 
corrections to the kinetic energy. 

The method can be extended to other fullerene systems of interest. It is a simple matter to 
treat endohedrally-doped fuUerenes or positively-charged fuUerenes. Lastly, by generalizing 
to finite temperature, equations of state can be calculated, as was done previously in the 
case of atoms [|TU|. 
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FIGURES 




FIG. 1. Electron charge density n (in units of ^ x 10^6 ^) vs a; for Z = 60 and R = 6.73ao. 
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FIG. 2. Total energy E vs shell radius R. 
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